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Abstract — Approximate  Dynamic  Inversion  (ADI)  has  been 
established  as  a  method  to  control  minimum-phase,  nonaffine-in- 
control  systems.  Previous  results  have  shown  that  for  single-input 
nonaffine-in-control  systems,  every  ADI  controller  admits  a  linear 
Proportional-Integral  (PI)  realization  that  is  largely  independent 
of  the  nonlinear  function  that  defines  the  system.  In  this  report, 
we  first  present  an  extension  of  the  ADI  method  for  single-input 
nonaffine-in-control  systems  that  renders  the  closed-loop  error 
dynamics  independent  of  the  reference  model  dynamics.  The 
equivalent  PI  controller  will  be  derived  and  both  of  these  results 
are  then  extended  to  multi-input  nonaffine-in-control  systems. 

Index  Terms — Dynamic  inversion,  feedback  linearization, 
proportion-integral  control,  PI  control. 

I.  Introduction 

YNAMIC  inversion  (DI)  or  feedback  linearization  is 
a  popular  control  design  method  that  is  well  suited 
for  minimum-phase  nonlinear  systems  [1]  [2,  Chapter  13]. 
DI  addresses  the  problem  of  designing  a  controller  to  trans¬ 
form  a  nonlinear  system  to  a  linear  one  by  feedback.  To 
overcome  some  limitations  imposed  by  the  requirements  of 
exact  linearization,  approximate  linearization  has  emerged  as 
a  viable  alternative,  where  the  problem  is  relaxed  to  enlarge  the 
class  of  admissible  controllers  [3],  A  notable  departure  from 
the  approximate  linearization  literature  is  [4],  where  tracking 
control  of  nonaffine-in-control  systems  are  considered. 

In  [4],  an  Approximate  Dynamic  Inversion  (ADI)  con¬ 
trol  law  was  proposed  that  drives  a  given  minimum-phase 
nonaffine-in-control  system  towards  a  chosen  stable  reference 
model.  The  control  signal  was  defined  as  a  solution  of  “fast” 
dynamics,  and  Tikhonov’s  Theorem  [2,  Theorem  11.2,  pp. 
439  -  440]  in  singular  perturbation  theory  was  used  to 
show  that  the  control  signal  approaches  the  exact  dynamic 
inversion  solution  and  that  the  system  states  approach  those 
of  the  reference  model  when  the  controller  dynamics  are  made 
sufficiently  fast. 


In  [5],  the  authors  showed  that  for  the  single -input  case, 
every  ADI  control  law  as  formulated  in  [4]  admits  a  linear 
Proportional-Integral  (PI)  model  reference  controller  realiza¬ 
tion.  The  key  characteristic  of  the  equivalent  PI  controller  is 
that  it  is  largely  independent  of  the  system’s  nonlinearities, 
in  contrast  to  the  original  ADI  control  law  in  [4],  However, 
when  the  controller  has  fast  dynamics  as  required  of  the  ADI 
method,  the  resulting  PI  controller  is  a  high-gain  controller 
with  associated  robustness  problems  [6]. 

Note  that  this  equivalence  holds  only  for  the  time  response 
when  applied  to  the  exact  system.  The  equivalence  do  not 
hold  when  applied  to  perturbed  systems  [7].  Even  when 
restricted  to  perfectly  known  minimum-phase  Linear  Time- 
Invariant  systems,  the  closed-loop  systems  differ  in  robustness 
properties  [7], 

This  report  extends  the  ADI  method  by  decoupling  the 
error  dynamics  specification  from  the  reference  model  dy¬ 
namics.  This  in  essence  decouples  the  “steady  state”  response 
specification  from  the  transient  response  specification,  when 
the  reference  model  response  is  viewed  as  the  “steady  state” 
response.  The  derivation  of  the  equivalent  PI  controller  for  this 
extension  is  very  similar  to  [5].  A  key  result  of  this  report  is 
to  extend  the  previous  statements  of  PI  and  ADI  equivalence 
to  multi-input  systems. 

The  report  will  proceed  as  follows.  Section  II  presents  the 
ADI  extension  and  PI  equivalent  controller  for  single-input 
systems.  The  second  (and  last)  section  extends  these  results  to 
multi-input  systems. 

II.  Equivalence  for  Single-Input  Systems 
A.  Nomenclature 

In  the  sequel,  italicized  symbols  (eg.  x)  denote  scalars, 
boldface  lowercase  letters  (eg.  x)  denote  column  vectors,  and 
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boldface  uppercase  letters  (eg.  A)  denote  matrices.  Upright 
text  subscripts  (eg.  xr  with  text  subscript  “r”  to  indicate  state 
of  reference  model)  are  variable  class  indicators,  and  italicized 
subscript  symbols  (eg.  xp  with  subscript  “p”  to  indicate  the  p- 
th  element  of  the  vector  x)  are  variables  for  numeric  quantities. 


B.  Approximate  Dynamic  Inversion  for  Single  Input  Systems 
Here,  the  ADI  method  [4]  for  single  input  systems  is  stated 
with  a  minor  generalization,  together  with  the  main  result.  The 
proof  in  [4]  applies  with  appropriate  (trivial)  substitutions,  and 
will  not  be  replicated  here. 

Consider  an  n-th  order  single-input  nonaffine-in-control 
system  of  relative  degree  p,  expressed  in  normal  form 


where 


x(f)  =  f(x(f),z(t),u(t)),  x(0)  =  x0, 

z(t)  =  g(x(t),z(t),u(t)),  z(0)  =  z0, 

x(f)  =  [x1(t),x2(t)>. . . ,  Xp(t)\T  G  Rp, 

r  x2(t)  i 


f(x(t),z(t),u(t))  = 


Xp(t) 

J(x(t),z(t),u(t))_ 


for  (x(t),z(t),u(t))  G  Dx  x  Dz  x  Du,  and  the  sets  Dx  C 
Rp,  Dz  C  M”_p  and  Du  C  ffi.  are  domains  containing  the 
origins.  Here,  [xT(t),  zT(f)]T  denotes  the  state  vector  of  the 
system,  u(t)  is  the  control  input,  and  f :  Dx  x  Dz  x  Du  i— >  R, 
g:  Ac  X  Dz  x  Du  i — >  R"  p  are  continuously  differentiable 
functions  of  their  arguments.  Furthermore,  assume  that  ^  is 
bounded  away  from  zero  for  (x(f), z(f), u(t))  G  fi  C  Dx  x 
Dz  x  Du,  where  is  a  compact  set.  That  is,  there  exists 
b0  >  0  such  that  ||£|  >  bo  for  all  (x(t),z(t),u(t))  G  fl. 
Note  that  |§£|  >  b0  >  0  implies  sign(|{)  G  {— 1,  +1}  is 
a  constant.  In  addition,  assume  that  the  function  /  cannot  be 
inverted  explicitly  with  respect  to  u. 

It  is  desired  for  x  (t)  to  track  the  states  of  a  stable  p-th  order 
linear  reference  model  described  in  the  controllable  canonical 
form 


xr(t)  =  Arxr(t)  +  Brr(f),  xr(0)  =  xr0,  (2a) 

where 


Here,  r(t)  is  a  continuously  differentiable  reference  input 
signal,  and  xr(f)  is  the  state  of  the  reference  model. 

Let  elf)  =  x(t)  —  xr(f)  G  Mp  be  the  tracking  error  signal, 
and  let  the  desired  stable  error  dynamics  be  specified  by 

e(i)  =  A  ce(t),  (3) 

where  Ac  is  Hurwitz  and  has  identical  structure  as  Ar,  but 
with  coefficients  aei  in  place  of  ar,  for  i  G  {0, 1, . . . ,  p  —  1}. 

Observe  that  in  [4],  Ac  was  set  equal  to  Ar,  while  in  the 
above,  an  independent  Hurwitz  matrix  Ac  can  be  specified. 
In  typical  applications,  Ar  and  Br  can  be  used  to  specify 
the  desired  system  response  to  excitation  r(t),  and  Ac  can  be 
used  to  independently  specify  the  desired  error  dynamics.  That 
is,  how  quickly  the  system  response  approaches  that  of  the 
reference  model.  Thus  the  preceding  is  a  slight  generalization 
of  the  ADI  as  formulated  in  [4], 

The  open  loop  (time-varying)  error  dynamics  are  then  given 
by  the  system 

e(t)  =  f(e(t) +xr(t),z(t),u(t))  -  Arxr(f)  -  Brr(f), 

(4) 

z  (t)  =  g(e(t)  +  xT(t),z(t),u(t)), 

with  initial  conditions  e(0)  =  eo,  z(0)  =  zq.  Define  the 
selector  vector 

c  =  [0,...,0,1]T  G  Rp. 

The  ideal  dynamic  inversion  control  is  then  found  by  solving 
for  u(t)  from  the  system  of  p  equations 

f(x(t),z(t),u(t))  -  Arxr(f)  -  Brr(f)  =  Ace(f). 

If  this  is  a  system  of  p  arbitrary  nonlinear  equations,  there  are 
in  general  no  solutions  of  u(t)  since  there  are  more  equations 
than  available  degrees  of  freedom.  However,  by  the  judicious 
choice  of  system,  reference  model  and  error  dynamics  repre¬ 
sentation,  the  first  (p—  1)  equations  are  automatically  satisfied. 
Therefore,  solving  the  above  p  equations  reduces  to  solving 
the  single  equation 

f(x(t),z(t),u(t))- cT(Arxr(f)+Brr(f))  =  cTAce(t),  (5) 


xr(f)  =  [£ri(f),£r2(f),  •  •  •  ,a;rp(f)]T  G  Rp,  (2b) 


and  the  Hurwitz  Ar  and  column  vector  Br  have  the  form 


resulting  in  the  exponentially  stable  closed-loop  tracking  error 
dynamics  (3).  Since  (5)  cannot  (in  general)  be  solved  explicitly 
for  u(t),  the  approximate  dynamic  inversion  controller  for  the 
above  formulation  can  be  given  in  similar  form  to  [4]  as 

eiiit)  =  —  sign  ^ f(t,e(t),z(t),u(t)),  (6a) 
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where 

f(t,e{t),z(t),u(t))  =  f{e(t)  +  xr(f),z(f),w(f)) 

(6b) 

-  cT(Arxr(f)  +  Brr(f)  +  Aee(f)), 

for  some  initial  control  it(0)  =  Uq.  Here,  e  is  a  posi¬ 
tive  controller  design  parameter,  chosen  sufficiently  small 
to  achieve  closed-loop  stability  and  approximate  dynamic 
inversion.  Observe  that  (6)  relaxes  the  requirement  for  exact 
dynamic  inversion  while  increasing  the  control  in  a  direction  to 
reduce  the  discrepancy  (5)  so  as  to  approach  the  exact  dynamic 
inversion  solution. 

Let  u  =  h(t,e,  z)  be  an  isolated  root  of  /(f,e,z,w)  =  0. 
In  accordance  with  the  theory  of  singular  perturbations  [2, 
Chapter  11],  the  reduced  system  for  (4)  is 

e(f)  =  Ace(t),  e(0)  =  e0, 

z(t)  =  g(e(£)  +  xr(i),  z(i),  h(t,  e(f),  z(f))),  z(0)  =  z0. 

With  v  =  u  —  h(t,e,  z),  and  r  =  t/e,  the  boundary  layer 
system  is 

^  =  -  sign  (j^j  f(t,e,z,v  +  h(t,e,z)).  (7) 

The  main  result  of  [4]  for  single-input  systems,  adapted  for 
the  generalization  above,  is  stated  below. 

Theorem  1.  Assume  that  the  following  conditions  hold  for  all 
(t,e,z,u  —  h(t,e,  z),e)  G  [0,  oo)  x  De  z  x  Dv  x  [0,eo]  for 
some  domains  De  z  C  and  Dv  C  R,  which  contain  the 

origins. 

1)  On  any  compact  subset  of  De  z  X  Dv,  the  functions  f 
and  g  and  their  first  partial  derivatives  with  respect 
to  (e,z ,u),  and  the  first  partial  derivative  of  f  with 
respect  to  t  are  continuous  and  bounded,  h(t,e,z) 
and  ^(t,e,z,u)  have  bounded  first  derivatives  with 
respect  to  their  arguments,  J^(£,  e,  z,  h(t,  e,  z))  and 
^(t,e,z,h(t,e,z))  are  Lipschitz  in  e  and  z,  uniformly 
in  t. 

2)  The  origin  is  an  exponentially  stable  equilibrium  of  the 
system 

z{t)  =  g (xr (f ) , z (f ) ,  h(t,  0, z (£))). 

The  mapping  (e,z)  i— >  g(e  +  xr(f),  z,  h(t,  e,  z))  is  con¬ 
tinuously  differentiable  and  Lipschitz  in  (e,z)  uniformly 
in  t. 

3)  ( t,e,z,v )  i— >  \^(t,e,z,v  +  h(t,e,z))\  is  bounded 
from  below  by  some  positive  number  for  all  (f,  e,  z)  G 
[0,  oo)  x  Dez. 


Then  the  origin  of  (7)  is  exponentially  stable.  Moreover,  let 
Q,v  be  a  compact  subset  of  Rv,  where  Rv  C  Dv  denotes  the 
region  of  attraction  of  the  autonomous  system 

dv  (  d  f  \ 

—  =  -  sign  f  —  J  /(0,eo,zo,u  +  /i(0,eo,Zo)). 

Then  for  each  compact  subset  fle.z  C  De  z,  there  exists  a 
positive  constant  e*  and  T  >  0  such  that  V  t  >  0,  (eo,zo)  G 
Oe,z)  u o  —  /i(0,eo,zo)  G  f lv,  and  V  e  G  (0,e»),  system  (1), 
(6)  has  a  unique  solution  xe(f)  on  [0,oo)  and 

x£(t)  =  xr(t)  +  O(e) 

holds  uniformly  for  t  G  [T,  oo). 

A  proof  of  Theorem  1  is  provided  in  [4],  In  summary. 
Theorem  1  states  that  when  regularity  assumptions  on  the 
system  dynamics  are  satisfied  to  ensure  existence  and  unique¬ 
ness  of  solutions,  and  system  (1)  is  minimum  phase  and 
controllable,  the  ADI  control  signal  u(t)  approaches  that  of 
the  exact  dynamic  inversion  solution,  and  the  system  states 
x(f)  approaches  and  maintains  within  ()(e)  of  the  reference 
model  states  xr(f)  for  a  sufficiently  small  e.  See  [4]  for  ways 
to  verify  the  assumptions  and  further  discussions. 

C.  Equivalent  PI  Controller 

Here,  we  recall  the  main  result  of  [5],  which  extends  trivially 
for  the  above  ADI  generalization.  For  notational  convenience 
in  the  sequel,  let  sign  ( )  =  a. 

Lemma  1.  For  every  Approximate  Dynamic  Inversion  con¬ 
troller  (6)  with  u{ 0)  =  Mo,  there  exists  a  linear  Proportional- 
Integral  model  reference  controller  realization 

u(t )  =  -^a(cTe(£)  -  g(t)),  (8a) 

where 

gif)  =  [  cTAce(A)  dX,  (8b) 

Jo 

and  g( 0)  =  cTe(0)  +  eauo- 

Proof:  Observe  that  from  (1),  we  have 

/(x(f),z(f),«(f))  =  ctx(£),  (9) 

and  from  (2),  we  have 

cT(Arxr(f)  +  Brr(f))  =  crxr(f).  (10) 

Substituting  (9)  and  (10)  in  the  ADI  control  law  (6)  gives 

eii(t)  =  —  acT(x(f)  —  xr(t)  —  A  ce(£)), 

=  — acT (e(£)  -  Aee(£)), 


(11) 
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which  can  be  integrated  in  time  to  yield  (8).  Setting  p(0)  = 
cTe(0)  +  ecrao  recovers  the  original  controller  initial  value 
u(0)  =  Uq.  ■ 

Note  that  a  =  sign(^)  G  {— 1,+1}  satisfy  a2  =  1. 
It  can  be  seen  that  the  result  (8)  is  a  PI  controller  acting 
on  the  error  between  the  system  states  and  the  states  of  the 
reference  model.  Furthermore,  observe  that  when  expressed  in 
the  error  coordinates,  e(£),  the  PI  controller  is  not  explicitly 
dependent  on  Ar  that  specifies  the  reference  model  dynamics, 
in  contrast  to  the  form  in  [5],  This  characteristic  is  the  result  of 
introducing  the  independent  matrix  Ac  for  the  error  dynamics 
specification.  Finally,  it  is  important  to  note  that  (9)  and  (10) 
are  true  by  definition,  while  (3)  is  a  design  specification 
that  the  controller  attempts  to  achieve,  so  that  (3)  cannot  be 
substituted  in  (11)  for  further  simplifications.  From  (8),  it  is 
apparent  that  the  PI  controller  attempts  to  achieve  (3),  which 
is  equivalent  to  achieving  (5). 

The  significance  of  Lemma  1  is  threefold: 

1)  The  PI  controller  allows  a  very  simple  exact  realization 
of  the  ADI  control  law. 

2)  The  PI  controller  is  a  linear  realization  of  a  (in  general) 
nonlinear  control  law. 

3)  The  PI  controller  realization  is  independent  of  the  non¬ 
linear  function  /(x(£),  z(£),  u(t))  in  (lb),  except  for  the 
sign  of  the  control  effectiveness,  a  =  sign(|£). 

The  existence  of  a  linear  realization  of  a  nonlinear  control  law 
hinges  critically  on  the  structure  of  the  underlying  system, 
reference  model,  error  dynamics  and  control  law.  Note  that 
the  PI  realization  does  not  apply  to  the  ADI  variant  in  [8] 
as  that  variant  uses  the  Jacobian  map  ^(x(£),z(£),w(£))  in 
place  of  sign(|£),  which  is  constant  by  assumption.  The  main 
practical  challenge  when  applying  the  method  of  [8]  lies  in 
obtaining  a  sufficiently  accurate  description  of  the  Jacobian 
map,  in  addition  to  the  nonlinear  function  /(x(£),z  (£),«(£)) 
in  (lb). 

As  an  aside,  adaptive  variants  to  [4]  have  been  proposed 
in  [9]— [  14]  in  which  the  nonlinear  function  /(x(£),z (£),«(£)) 
is  assumed  unknown.  These  variants  attempt  to  estimate 
the  unknown  function  /(x(£),  z(£),  «(£)),  and  construct  an 
analogous  ADI  control  law  based  on  the  estimate.  Since  the 
PI  controller  realizes  the  ADI  control  law  exactly  without 
explicit  dependence  on  /(x(£), z(£), it(£)),  it  appears  that 
these  approximation  based  adaptive  variants  are  unnecessary. 
Numerical  results  in  [5],  [15]  show  that  for  the  single¬ 
input  case,  the  PI  controller  achieves/exceeds  the  tracking 
performance  of  the  adaptive  variants  proposed  in  [13],  [14] 


respectively. 


III.  Equivalence  for  Multi-Input  Systems 

Here  we  extend  the  above  results  to  multi-input  nonaffine- 
in-control  systems. 


A.  Approximate  Dynamic  Inversion  for  Multi-Input  Systems 


Here,  the  ADI  method  [4]  for  multi-input  systems  is  stated 
with  a  minor  generalization  analogous  to  the  single  input  case. 
Consider  a  multi-input  nonlinear  system  expressed  in  normal 
form 


x(£)  =  f(x(£),z(£),u(£)),  x(0)=x0, 

z(£)  =  g(x(£),z(£),u(£)),  z(0)  =  z0, 

where 

u(£)  =  [u1{t),u2{t),...,um{t)]r  £  Rm, 
x(£)  =  [x?(£),xj (£),..., x^(£)]T  g  Rp, 

fi(x(£),z(£),u(£)) 

,  ,  ,  .  ,  ...  f2(x(£),z(£),u(£)) 

f(x(£),  z(£),  u(£))  = 

fm(x(£),z(£),u(£))J 

and  for  k  £  {1, 2, ... ,  m},  with  Pk  =  p, 

xfc(£)  =  [xki(t),xk2(t),  •  •  •  ,XkPk{t)\T  G  1 

Xk2{t) 


ffc(x(£),  z(£),  u(£))  = 


Xkpk  (£) 

L/fc(x(f)>  z(£),  u(£))J 


(12a) 


(12b) 


(12c) 


Here,  [xT(£), z1  (£)]T  G  R"  is  the  state  of  the  n-th  order 
system  with  m  inputs,  and  u(£)  is  the  vector  of  control  inputs. 

The  objective  is  to  design  u(£)  so  that  the  state  x(£)  tracks 
the  state  xr(£)  of  the  stable  p-th  order  reference  model  in 
Brunovsky  canonical  form 


xr(£)  =  Arxr(t)  +  Brr  (£),  xr(0)  =  xr0, 

where 


r(£)  =  br(£),r2(£),  •  •  •  ,rm{t))T  G  Rm, 
xr(£)  =  [xJ1(£),x^2(£),...,xJm(£)]T  GRp, 


Ari 

Bri 

,  B,  = 

A 

■^rm 

^rm 

(13a) 


(13b) 
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and  for  k  G  {1,2,...,  to},  with  pk  as  in  (12), 
Xrfc(f)  =  [ Xrkl{t) ,  Xlk2{t ),  •  •  •  ,Xrkpk(t)]T  G  RPfc, 


0 

l  •• 

0 

'  0  ' 

0 

o  •• 

1 

1  B pfc  - 

0 

_  ^rkl 

—Ork2  '  ' 

~arkpk_ 

J^rk_ 

(13c) 


Here,  each  Arfc  for  k  G  {1,2,...,  to}  is  a  specified  Hurwitz 
matrix,  r(t)  is  the  vector  of  continuously  differentiable  refer¬ 
ence  inputs,  and  xr(£)  is  the  state  of  the  reference  model. 


The  tracking  error  vector  can  be  defined  and  decomposed 
in  correspondence  with  the  above  system  and  reference  model 
state  variables  as 

e(t)  =  x(t)  —  xr(t) 

=  [e? (t),ej(t ),..., e^(f)]T  G  Kp, 
where  for  he  {1,2,... ,  to}, 

ek(t)  =  xfc(f)  -  xrfc(f) 

=  [ekl (t),ek2(t), . . .  ,ekpk(t)]T  G  KPfc, 
and  for  each  k  G  {1,2,...,  to},  i  G  {1,2,...,  pk}. 


or  equivalently,  solving  the  system  of  to  scalar  equations 

/fc(x(t),z(t),u(t))  -cj (A rfcxrfe(f)  +Brfcrfc(f)) 

t  (16) 

for  k  G  {1,2,..., to},  resulting  in  the  exponentially  stable 
closed-loop  tracking  error  dynamics  (14).  Since  (16)  cannot 
(in  general)  be  solved  explicitly  for  u(f),  the  approximate 
dynamic  inversion  controller  for  the  above  formulation  can 
be  given  in  similar  form  to  [4]  as 

eu(f)  =  Pf(f,e(f),z(f),u(f))  (17a) 

where  P  G  Rmxm  is  a  controller  parameter,  and  for  k  G 
{1,2, . . .  ,m},  the  fc-th  element  of  f(f,e(t),z(t),u(f))  is 
h(t,e(t),z(t),u(t))  =  fk(e(t)  +  xr(*),z(J),u(t)) 

(Ar^Xr/c(t)  +  +  A ek&k(t)')  i 

for  some  initial  control  u(0)  =  Uq.  Observe  that  for  the  multi¬ 
input  case,  there  is  no  analogue  to  the  sign  of  the  control 
effectiveness,  sign  (pp),  and  the  controller  design  parameters 
are  e,  P,  Ark,  Ack  and  Br(.  for  fc  G  {1,2,...,  m}. 

Let  u  =  h(f,e,z)  be  an  isolated  root  of  f(t,  e,z,u)  =  0. 
The  reduced  system  for  (15)  is 

e(t)  =  A0e(f),  e(0)  =  eo, 


Let  the  desired  stable  error  dynamics  be  specified  by 


z(t)  =  g(e(£)  +  xr(f),  z(f),  h(f,  e(f),  z(f))),  z(0)  =  z0. 


e(f)  =  Ace(f),  (14) 

where  Ac  has  identical  structure  as  Ar,  but  with  block  entries 
Ack  in  place  of  Alk  for  k  G  {1,2,...,  to},  each  Aok  is  chosen 
to  be  Hurwitz,  and  has  identical  structure  as  Ark,  but  with 
coefficients  aeki  in  place  of  arki  for  i  G  {1,2,...,  pk}.  Similar 
to  the  single-input  case,  this  extension  allows  the  desired  error 
dynamics  to  be  specified  independently  of  the  reference  model 
dynamics. 

The  open  loop  (time-varying)  error  dynamics  are  then  given 
by  the  system 

e(t)  =  f(e(i)  +  xr(f),z(f),u(f))  -  Arxr(£)  -  Brr(f), 
z(f)  =  g(e(f)  +xr(f),z(f),u(f)), 

with  initial  conditions  e(0)  =  eo,  z(0)  =  z(J.  For  k  G 
{1,2,...,  to},  let 

ct  =  [0,...,0,l]TGlp‘. 

The  ideal  dynamic  inversion  control  is  found  by  solving  for 
u  (t)  from  the  system  of  p  equations 


With  v  =  u  —  h(f,e,z),  and  r  =  t/e,  the  boundary  layer 
system  is 

^=Pf(t,e,z,v  +  h(t,e,z)).  (18) 

dr 

The  main  result  of  [4]  for  multi-input  systems,  adapted  for  the 
above  generalization,  is  stated  below. 

Theorem  2.  Let  the  following  conditions  hold  for  all  (t.  e,  u  — 
h(f,  e,  z),  e)  G  [0,  oo)  x  De  z  x  Dv  x  [0,  eo]  for  some  domains 
De  z  C  K"  and  Dv  C  Km  that  contain  the  origins. 

1)  On  any  compact  subset  of  De  z  x  Dv,  the  functions 

f  and  g  and  their  first  partial  derivatives  with  respect 
to  (e,z,u),  and  the  first  partial  derivative  of  f  with 
respect  to  t  are  continuous  and  bounded,  h(f,  e,z) 
and  z,u)  have  bounded  first  derivatives  with 

respect  to  their  arguments,  §§(t,  e,  z,  h(£,  e,  z))  and 
^  (t,  e ,  z ,  h (t,  e,  z))  are  Lipschitz  in  e  and  z,  uniformly 
in  t. 

2)  The  origin  is  an  exponentially  stable  equilibrium  of  the 
system 


f(x(f),z(f),u(f))  -  A rxr(f)  -  Brr(f)  =  Aee(f), 


z(£)  =  g(xr(f),z(f),h(f,0,z(f))). 
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The  mapping  (e,z)  i— >  g(e  +  xr(f), z,  h(f,  e, z))  is  con¬ 
tinuously  differentiable  and  Lipschitz  in  (e,z)  uniformly 
in  t. 

3)  For  every  (t,  e,z)  £  [0,  oo)  x  De  z,  all  the  eigenvalues 

of 


and  from  (13),  we  have 

cj  (A  rfcxrfe(f)  +  Brfcrfc(i))  =  c  £xrfc(f).  (21) 

Substituting  (20)  and  (21)  in  the  multi -input  ADI  control 
law  (17)  gives 


P 


dt , 

at(,-e'z-v 


h(£,e,z)) 


/juve  negative  real  parts  bounded  away  from  0. 


Then  the  origin  of  (18)  is  exponentially  stable.  Moreover,  let 
fiv  be  a  compact  subset  of  Rv,  where  Rv  C  Dv  denotes  the 
region  of  attraction  of  the  autonomous  system 

<T~  =  Pf(0,eOlzo,v  +  h(0,eo,zo)). 
dr 

Then  for  each  compact  subset  fle,z  C  De  z,  there  exist  a 
positive  constant  e*  and  T  >  0  such  that  Vi  >  0,  (eo,zo)  £ 
fie,z>  uo  —  h(0,eo,zo)  £  fiv>  and  Ve  €  (0,  e*),  system  (12), 
(17)  /zas  a  unique  solution  xe(i)  on  [0,  oo)  and 


xe(f)  =  xr(t)  +  O(e) 


/zoMs  uniformly  for  t  £  [T,  oo). 

In  essence.  Theorem  2  is  the  multi-input  extension  of 
Theorem  1.  See  [4]  for  ways  to  verify  the  assumptions  and 
further  discussions. 


B.  Equivalent  PI  Controller 


eii(f)  =  P 


c?  (xi(f)  -  xri(f)  -  Aciei(f)) 
cj (x2(f)  -  xr2(i)  -  Ac2e2(i)) 


|_cm(xm(£)  xrm(t)  Aemem(£))  J 

"  c?(®iW  “  Aeiei(i)) 

(e2(i)  -  A c2e2(i)) 

.  ? 

Aem^m(i))_ 


which  can  be  integrated  in  time  to  yield  (19).  Setting  w(0)  = 
v(0)  —  <=P 1  U(,  recovers  the  original  controller  initial  value 
u(0)  =  u0.  Since  P  must  be  chosen  to  satisfy  assumption  3 
of  Theorem  2,  all  eigenvalues  of  P  [  J^-]  have  strictly  negative 
real  parts,  which  implies  that  P[Jjj]  have  full  rank.  This  in 
turn  implies  that  P  has  full  rank  which  ensures  the  existence 
ofP-1.  ■ 

In  (19),  the  vector  v(£)  is  formed  from  elements  of  the  error 
vector  e(i),  and  w(£)  is  formed  from  linear  combinations  of 
the  integral  of  the  error  vector.  It  can  be  seen  that  (19)  is  a 
multi-input  PI  controller  that  realizes  the  ADI  controller  (17). 
Similar  to  the  single -input  case,  we  note  that  the  PI  realization 
does  not  apply  to  the  ADI  variant  in  [8]. 


Here,  we  extend  the  previous  results  to  show  that  for  the 
multi-input  case,  there  also  exists  an  equivalent  PI  controller. 

Lemma  2.  For  every  Approximate  Dynamic  Inversion  con¬ 
troller  (17)  with  u(0)  =  Uo,  there  exist  a  linear  Proportional- 
Integral  model  reference  controller  realization 


u(f)  =  ~P  (v(f)  —  w(f)) 
e 

where 

v(f)  =  [cj'ei (t ), c Je2 (t ),..., c^em(f)]T  € 

r  So  C1  Aelel(/^)  d\ 

fo  c2  Ae2e2(A)  dX 


r(t)  = 


L/0  VAj 


(19a) 


(19b) 


and  w(0)  =  v(0)  —  eP  1u0. 

Proof:  For  k  £  {1,2,...,  m},  we  have  from  (12), 


/fc(x(f),z(f),u(f))  =C  feXfc(f), 


(20) 


IV.  Conclusions 

An  extension  of  the  Approximate  Dynamic  Inversion  (ADI) 
method  for  minimum-phase  nonaffine-in-control  systems  was 
presented  that  renders  the  error  dynamics  independent  of 
the  reference  model  dynamics.  In  essence,  this  decouples 
the  “steady  state”  response  specification  from  the  transient 
response  specification,  where  the  “steady  state”  response  is 
specified  by  the  reference  model  dynamics  while  the  transient 
response  is  independently  specified  by  the  error  dynamics.  It 
was  shown  that  every  ADI  control  law  admits  an  equivalent 
linear  Proportional-Integral  (PI)  controller  realization  that  is 
largely  independent  of  the  nonlinearities  of  the  system,  for 
both  single-input  and  multi-input  systems. 

We  note  that  this  equivalence  must  be  interpreted  with 
caution.  As  shown  in  [7],  this  equivalence  holds  only  for  the 
time  response,  and  only  when  applied  to  the  exact  system.  In 
particular,  even  when  specializing  to  minimum-phase  linear 
time-invariant  systems,  they  differ  in  robustness  properties. 
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